Radiation Due to Josephson Oscillations in Layered Superconductors 
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We derive the power of direct radiation into free space induced by Josephson oscillations in intrinsic Joseph- 
son junctions of highly anisotropic layered superconductors . We consider the super-radiation regime for a crystal 
cut in the form of a thin slice parallel to the c-axis. We find that the radiation correction to the current-voltage 
characteristic in this regime depends only on crystal shape. We show that at large enough number of junctions 
oscillations are synchronized providing high radiation power and efficiency in the THz frequency range. We 
discuss crystal parameters and bias current optimal for radiation power and crystal cooling. 
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Josephson junctions (JJs), as sources of tunable continuous 
electromagnetic radiation, were discussed for a long time after 
the prediction of the ac Josephson effect.|l] Early measure- 
ments 1 2] demonstrated that emittence from a single J J has 
very low power, typically ^ 10^*' /iW. Since then a signifi- 
cant effort has been devoted to develop JJ arrays as coherent 
sources of radiation, see, e.g., Refs. Q, 0. A major challenge 
is to force all JJs in array to emit coherently, so that power 
increases proportionally to the square of the total number of 
junctions. (sll In particular, for an array of 500 junctions a max- 
imum power of the order of 10 /iW at discrete frequencies 
< 0.4 THz has been achieved so far in the super-radiation 
regime. [5] The difficulties to synchronize many artificial JJs 
are related mainly to the facts that artificial junctions always 
have slightly different parameters, especially the Josephson 
critical current, and that one cannot put many of them at dis- 
tances smaller than a wavelength but needs to distribute them 
over a wavelength or more.yl \§\ Also, as the maximum fre- 
quency is limited by the superconducting gap, it can not ex- 
ceed several hundred gigahertz for structures fabricated out of 
conventional superconductors. 

Layered high-temperature superconductors like 
Bi2Sr2CaCu208 (BSCCO) offer a very attractive alter- 
native for developing radiation sources. ja S S S A large 
value of the gap (up to 60 meV) allows for very high Joseph- 
son frequencies, which can be brought into the practically 
important terahertz range. Moreover, intrinsic JJs (IJJS) 
have much closer parameters than artificial ones as these 
parameters are controlled by the atomic crystal structure 
rather than by amorphous dielectric layer in artificial JJs. 
Also, layered superconductors provide a very high density 
of IJJs (1 per 15.6 A along the c-axis) and thus it is easy 
to reach super-radiation regime with many junctions on the 
scale of radiation wave length. In this regime the radiated 
electromagnetic field effectively couples JJs and helps to 
synchronize them. In this Letter we demonstrate that the 
superradiation regime indeed results in the synchronization 
of Josephson oscillations in IJJs in zero external magnetic 
field. Thus c-axis current biased crystal may work as source 
of Josephson coherent emission of radiation (JOCER). We 
calculate the radiation power and IV characteristics and 




FIG. 1 : (color online) Left: Schematic picture of layered supercon- 
ductor placed in between plate-like leads serving as screens. The 
directions of the dc transport current J, and of the radiation Poynt- 
ing vectors Px are shown. Right: Influence of radiation on current- 
descending branch of the IV dependence. 



discuss an optimal crystal geometry accounting for heating 
due to quasiparticle dissipation. 

So far mostly radiation from the flux flow of the Josephson 
vortices was discussed in the literature. The induc- 
tive interlayer coupling typically promotes formation of the 
triangular vortex lattice. However, to generate noticeable out- 
side radiation oscillations induced by the moving lattice have 
to be in phase in different layers, which is realized only if the 
moving vortices form a rectangular lattice. No regular way 
to prepare such a lattice is known at present. In addition, it 
seems to be unstable in most of parameter space. HEl Here 
we consider the synchronization of the Josephson oscillations 
by radiation field in the simplest case, when dc magnetic field 
is not applied and only radiation itself introduces the in-plane 
phase gradients. 

In the resistive state phases ipn oscillate at the Josephson 
frequency, ujj = 2eV/h, where the voltage V between the 
neighboring layers is induced by interlayer dc current. For 
uniform oscillations and identical junctions the voltage V is 
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the same in all junctions (except, possibly, top and bottom 
junctions in the stack) because the same current flows between 
all layers. To reach a resistive state the pulse of dc current ex- 
ceeding the Josephson critical current should be applied but 
then current may be diminished to reach necessary voltage. 
Transport measurements in BSCCO have shown that the resis- 
tive state on current-descending branch is preserved down to 
the voltage Vp ^ hu!„ /2e, where Wp, is the Josephson plasma 



frequency see, e.g., Ref . 1 1 3L [u!p/(2Tr) « 0.15 THz in opti- 
mally doped BSCCO]. We consider crystal with sizes Lx, w 
and Lz ~ Ns in the directions x,y, and z respectively, see 
Fig. [U The conditions Lx,Lz < — c/ujj, are neces- 
sary for the super-radiation regime. The boundary conditions 
for the oscillating phase are sensitive to distribution of outside 
em fields which, in turn, depend on geometry of the stack and 
electric contacts. We consider the simplest geometry assum- 
ing that (i) w ^ k^^ so that all quantities are y independent, 
and (ii) JJ stack is bounded by metallic or superconducting 
contacts with the same lateral sizes as the stack which extend 
in the z dkection over distance Lgc ^ fcj^ (see Fig.[T]). We 
also assume that the contact material has a very small surface 
impedance so that the ac electric field at the contact surface is 
negligible. Such contacts serve as "screens", restricting radia- 
tion to half-infinite spaces \x\ > Lx/2. This greatly simplifies 
the analytical consideration. As dissipation increases with 
and radiation is not, we argue that optimal crystal should be 
platelike with < L^- 

To find the phase differences ifin{x, t) inside the crystal we 
solve finite-difference differential equations, see, e.g., Ref.ll 1 



d^ifn / 2 \ f 9ipn ■ 
-g:^ = ("^n- 1) I Sinipn- Vuhy^n ) , (1) 



(2) 



- r'Tab hy,n + Tab V u^n = 0. 



We use reduced x coordinate, u ~ x/Xj normalized to 
= IS, reduced time, r = ujpt and u = uij/uip, with 
^p = c/ (Ac Y^), and reduced magnetic field hy^n = By^n/Bc 
with Be — <i>o/(27rAafcAc), where By^n is the magnetic field 
between the layers n and n+1. Here ec is the c-axis dielec- 
tric constant inside the superconductor, Xab and Ac are the 
London penetration lengths, 7 = \c/\ab is the anisotropy 
ratio. In terms of these parameters the Josephson critical cur- 
rent is Jc = $oc/(87r^sAc). Further, fab = l + Vabd/dr, 
I = Xab I s , V^j notates the discrete second derivative opera- 
tor, V^^An — An+i +An-i — 2 Aj, and a ~ 0. 1—1 is the parame- 
ter of the capacitive coupling. 11211 The dissipation parameters, 
Vab — 'i'^o-ab / {I'^^c'-^p) ^nd Vc = 47r (Tc/ (ccLdp), are determined 
by the quasiparticle conductivities, aab and ctc, along and per- 
pendicular to the layers, respectively. The electric field inside 
the superconductor between the layers n and n + 1 is given by 



(1 - a\7l)E,n = {BJ/^){d^n/dT). 



(3) 



The typical parameters of optimally doped BSCCO at low 
temperatures are ec = 12, s = 15.6 A, j — 500, Aab = 200 nm, 
Jc = 1700 A/cm2, crc(O) = 2 • lO^^ (ohm-cm)-\ aab{0) = 



4 • 10'' (ohm-cm)-i.|[ll This gives £ « 130, i^ab ~ 0.2, and 
i^c ~ 2 • 10~^. An important feature of BSCCO is higher 
in-plane dissipation in comparison with c axis one. 

The boundary conditions, i.e., relations between time and 
space derivatives of (p„ at the edges parallel to (y, z), are 
determined by the relations between the electric and mag- 
netic fields in the outside media. As the y and z sizes of 
the system (crystal and screens) are assumed to be larger than 
the wavelength, the dielectric media can be treated as infi- 
nite in these directions. Such a half-infinite space geom- 
etry allows us to find the boundary conditions analytically. 
From the Maxwell equations in the free space we find rela- 
tion between the magnetic, B = (0, By, 0), and the electric, 
E — {Ex,{),Ez), fields at the boundaries. We assume that 
there are only outgoing waves from the crystal [y, z) edges 
meaning that the fields have the coordinate and time depen- 
dence exp(ifca;|a;|-HfczZ— iwr), where fc^; =sign(a;)(A;5— ^z)^^^ 
for k1 < k'^ and k^ = i{kl—k'^y/^ for kl>k'^. The relations 
between fields at u~±Lx/2 (Lx = L^/ Xj) are ifTsIl 



By{uj, kz) 
CJkz) = 



TCuj{kz)Ez{Lu,kz), 
\k^\{kl-kl)-^/\ for kl<kl, 
-ik^{kl - kl)-^l^, for kl > k^. 



(4) 



Inverse Fourier transform with respect to k^ gives nonlo- 
cal relation between the magnetic and electric fields at the 
edges. As we assume that the screen material has small sur- 
face impedance, we can neglect the electric field at |z| > Lz/2 
and, using Eq. (O, we obtain the reduced boundary condi- 
tion connecting the magnetic field /iy „ with the phases at the 
edges {hy,n is expressed via V(^„ by Eq. (|2])), 

isiuj v-^, 2 ^-l 

X [\k^\Jo{k^s\n~m\) + ik^NQ{k^s\n-m,\)], (5) 

where Jo (a;) and No{x) are the Bessel functions. 

We consider high-frequency Josephson oscillations, uj ~ 
ujj/ujp ^ 1, in the case of layered crystals with large num- 
ber of junctions N > £. This allows us to neglect finite-size 
effects along the z-axis. The equation for uniform solution 

ip„{u,T) = ip{u,T) is 



aV/Sr^ + vrdif/dr + sin - f'^lip = 0. 



(6) 



In the limit w 3> 1 we look for the solution in the form 
</j(u, r) = ujT + (f){u, t) with 0^1. Eq. (|5]l gives the bound- 
ary conditions for (paiu — ±.Lx/2 ifTsIl 



(7) 
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The solution is r) = Im [(t)^{u) exp(— iwr)], where 

(l)^ — — {ijj'^ + iijjVc)^^+Acos{kuju) \N\\hkuj—ijj/£, (8) 
A^iC,{^ui s\-n(kujLx/2) + iC,u)COs{ku,Lx/2)]{uj + iVc)}^^ ■ 
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Here \(\ <C 1. The first term in (j)^^ is the amplitude of Joseph- 
son oscillations, while the second term describes the electro- 
magnetic waves propagating inside the junctions. They are 
generated at the boundaries due to the radiation field. 

Next we show that coherent radiation field, similar to all 
junctions, in combination with intralayer dissipation stabilizes 
the uniform Josephson oscillations. For that we have to con- 
sider a small perturbation to the uniform solution, (pn{u, t) = 
ujT+(j){u, r r) and verify that there is no perturbations, 

'i9„(u, t), increasing with time. The general analysis is rather 
cumbersome. Results in a closed form may be obtained only 
by use of approximations valid in the limiting case considered 
here, 1^,1^ < fcj\ > 1, and ^ > 1. 

Equations for {u, r) are obtained by linearization of Eqs. 
^ and (|2]i. The term cos[(/i(t)]i?„(m, r) « cos(tJT)z9„(M, t) 
in the linearized equation couples harmonics with small fre- 
quency n with the high-frequency terms ft ± u. At uj ^ 1 
we can neglect coupling to the higher frequency harmonics 

± mu! with m > 1 and represent the phase perturbation 
(and field) as 



' E exp(i/3wr) 

/3=±1 



sin(gn) exp (^iflr) 



with q = T:k/{N + 1), k = 1,2 . . . , N. Here the complex 
eigenfrequency O = 0(g) is assumed to be small, <C 
and has to be found. Stability means that Im[il] < for all q. 
Substituting this presentation into the linearized equations ([U 
and (|2]i, excluding oscillating magnetic fields, and separating 
the fast and slow parts, we obtain coupled equations 



^2 



l + an 



(9) 
(10) 



Here C = (cos < 



« aq = aq^ with q^ = 2(1- 

cosq), Gl0 = qy[l-iin- f3u;)i^ab]+r^ and - G,.o. 
Using Eqs. ©, Q (in the Hmit fc^ > due to fc^ix < 
tt), and (|5]l we get the boundary conditions for slow and fast 
components at u = ±Lx/2 for q^ tt/N, 



(l + a,)ec97 



(11) 
(12) 



Because of the condition jOj <C t^-", in most cases one can ne- 
glect Vl in equation and boundary conditions for dq^p. We also 
assume Vc 1 ^ oj and neglect dissipation when it is not es- 
sential. As dq varies at the typical length scale ~ l/G^O, 
which is much larger than L^, the coordinate-dependent part 
of dq can be treated as a small perturbation. Neglecting the 
coordinate dependence of dq in the equation for d, we obtain 
the approximate solution of Eqs. ( fTOl i and ( fT2b . Substituting it 



into Eq. ([9]), we obtain Mathieu equation for the slow-varying 
component 

02 



l + aa 



+ t^^n-^-V{u) + G-'^i]dq^O, (13) 



where the "potential" is given by V{u) = Vi (u) + V2 {q, u), 



Vi(u) 



V2{q, — Re 



1 

2^ 
1 

2w 



Re 



k^ sm{k^Lx/2) + i(^uj cos{k^Lx/2) 
Kjf- cosp^u 



p+ sin(p+L^/2) + K+ cos(p+La;/2) 



and = wGg_+. In the super-radiation regime, kuiL^ — 
LoLx/t ^ 1, in the lowest order with respect to ix}Lx/(-, the 
part Vi{u) reduces to a constant, Vi{u) « K,^l(2uP') with 
= [/;^(/:^+eca) + l)/[(£„+eca)^ + l] sxyAa = LxlL^. 
Equation ( flj] ) and the boundary conditions (fTTI) determine 
the spectrum of small perturbations to the uniform solution. 
Treating the coordinate-dependent part of dq as a small per- 
turbation allows us to derive the expression for 0(5), 



W2{q) =Re [2/[p+i,(p+/K++cot(p+i,/2))] 



(14) 



From this result we can conclude that the main contribution to 
stabilization of uniform oscillations comes from the term K.^, 
describing effective coupling of junctions due to the radiation. 
Its stabilization effect increases with as K.^^ ~ C^^Lz/ tcLx 
for Lz < icLx and /C^ 1 for > £cLx- The charging- 
effect term, aq, also contributes to stabilization. The term W2 
describes the effect of modes dq,± induced inside the crystal 
due to radiation. Formally, the W2 term leads to instabilities 
in the limit of zero dissipation because its denominator van- 
ishes near the resonance values of q given by 2(1 — cos g) — 
[2Trm/{ujLx)]'^ ^ l/i"^, where m is an integer. These insta- 
bilities correspond to parametric excitation of the Fiske reso- 
nances described by Eq. ( fTOl ). However, they are suppressed 
already by very small dissipation. Indeed, at small dissipation 
we estimate the maximum value of W2 as [q^Sc lui[p^]Lx]^^, 
where Im[p+] « sin(g/2) (I'c + i'abi^^). As I'abW^ ^ I'c, we 
see that |11^2,max| < 1 for i/ah > l/(27r27ec) ~ 10"^. For real- 
istic level of dissipation in BSCCO, Vab ^0.2, the resonance 
features in W2 are completely washed out and IW2I ^ 1 for 
all g's. Thus the intralayer dissipation stabilizes uniform os- 
cillations but it does not affect them in any other way. 

We proceed now with derivation the radiation power and 
IV characteristics in the super-radiation regime of uniform 
Josephson oscillation. The Poynting vector at x = ±Lx/2 
in terms of the oscillating phase is given by jlSll 



PxiL0)^± 



(ft2, ,3 



jQ{k^s\n-'m\)\Lp^{±Lxl2)'^ , 



where the oscillating phase difference, ip,^, is determined by 
Eq. (O. Substituting this solution in the limit k^Lx ^ 1, 



4 



we obtain for the total radiation power Vra.d (w) = Px {lli)LzW 
going from one side. 



Cia)^a\l/[{ae, + C^f + 1], 



(15) 



For small L^, ^ Lz/^c, Vrad « L'^ and it is A^- 
independent, while for larger the geometrical factor 
C{a) 1 , meaning that T^rad JV"^ and it is iv-jj-independent. 

The dc interlayer current density J consists of quasiparticle 
contribution, acV/s, and the Josephson part Jc sin ip, aver- 
aged over time and coordinate. We derive for j — J/ Jc, 



J = VcUJ + - Im 0^ )„ = ly^uj + —3 + — ^ . 



(16) 



The last term, jiad, describes contribution to the dc current 
due to radiation losses. This part of the current multiplied 
by the total voltage gives total radiation power 27-'iad- As 
a function of uj the current j has minimum at = ujai ~ 
[C{a)/eci'ca)]^^^ ^ 1. The influence of radiation on the IV 
dependence is illustrated in Fig.[T] Only part of the IV char- 
acteristics at (jj > w„i is stable. At the voltage V„j — fiw„i/2e, 
corresponding to the current j„i — (3/2)z^c^m^ the stack jumps 
back to the static state. At this "retrapping" current the dissi- 
pation power is twice of the radiation power, i.e., for the con- 
version efficiency we obtain 2Vrad/ {Vdis + SPiad) < 1/3. 

An important issue is the stability of the coherent state with 
respect to parameter variations from layer to layer, which may 
include the crystal width, 6Lxn, and the Josephson current 
density, 6 Jen- The coherent solution with the same voltage 
drop in all junctions, V, must satisfy the current conserva- 
tion condition. Such solution can be built as ipn{u,T) = 
i^T + (f>^^n{u,T) + Pn, wherc the additional phase shifts, /3„, 
compensate for parameter variations. In the case of smooth 
parameter variations and small charging parameter, we can 
derive correction to the local current and obtain 



sin(/3„) - S : 



SL 

xn 



{Lxn) Jrad {Jen) 



with S = (sin(/3„)) and (...) means average over n. We 
can see that the coherent state survives until SLxn/ (Lxn) < 
jia.d/{vc^^) and 5 Jen < {Jen)- These conditions do not im- 
pose too demanding restrictions on the acceptable range of 
parameter variations. 

Evaluating the Joule heating, we find that the cooling rate 
Q per unit area of each crystal side || yz should be 



Q 



327r3A2s2 



Lx 



to 



27^rad 



(17) 



The maximum efficiency is reached at a; ~ lo^. To achieve 
this at frequencies close to 1 THz one needs to maximize 



cx [Cia) /a] "'^/^ by optimizing the crystal shape. At ~ 1 
we get max[£(a)/a] « ec/4 at a w I/ec- This gives ~ 5 
corresponding to ijJj/{2ti) = 0.75THz. At this frequency, as- 
suming Lz — 40 /im, the optimum lateral sizes are L^; « 4 /im 
and w > 300 /im. Biased with the current density « O.OlJc, 
such a crystal radiates with the power T'rad/w ~ 30mW/cm 
from each side, while it should be cooled with the rate Q w 15 
W/cm^ at each side. 1 16] As Lx increases, Vi-ad/w saturates 

— 1/3 

to w 0.2 W/cm at u! Um Lx , while Q increases lin- 
early with Lx - We note that increasing the number of layers 
also promotes synchronization of oscillations by the radiation 
field. 

In conclusion, we have shown that uniform Josephson os- 
cillations in intrinsic junctions of layered superconductors 
are stable in the superradiation regime at high frequencies 
ujj ^ ujp. They lead to coherent radiation into free space with 
significant power and efficiency as high as 1/3 at frequencies 
^ 1 THz. Important point is that to have reasonable cooling 
rate and strong radiation the crystal should be in the form of 
thin plate along the c with large number of layers N ^ 10''. 
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